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INTRODUCTION 




Over the past three decades, much effort has been devoted to solving the 
inviscid flow over blunt bodies. Two basic approaches have been used: (1) 

an inverse approach - where the shock wave is given and the body shape and 
details of the flow within the shock layer are unknown and (2) a direct 
approach - where the body shape is given and all details of the flow within 
the shock layer are unknown. (See ref. 1.) 

i 

For the inverse problem, the steady flow equations are usually solved by 
marching inward from the shock wave to obtain the solution for the body. 

(See refs. 2 to 5.) Since the steady flow equations are of a mixed 
elliptic-hyperbolic type, this initial value problem is ill posed and is 
inherently unstable. However, for some cases good results can be obtained 
by this approach before the instabilities destroy the solution (refs. 1 and 
2). The inverse problem can be made direct by iteration. However, 
improvements to a solution which gives an approximation to a desired body 
shape is not always apparent (ref.1), and thus solutions for complex body 
shapes are usually difficult to obtain. 

More recently, time-dependent, finite-difference techniques have been 
used with great success to solve the more difficult direct problem. (See 
refs. 6 to 11). This technique takes advantage of the hyperbolic nature of 
the unsteady flow equations by marching in time from an approximate initial 
condition to the desired steady-state solution. This problem is well posed 
and thus eliminates the instabilities that were mentioned earlier for the 
inverse problem. This technique has been used successfully to compute the 
flow over very complex geometries. 


The purpose of the present paper is to describe such a time-dependent, 
finite-difference computer code, developed by the senior author, which has 
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fields over axisymmetric bodies at zero degrees angle of attack. The code 
has the capability of calculating the flow of an ideal gas, CF^, or air in 


chemical equilibrium and can be applied to reasonably general axisymmetric 
body shapes. The present paper gives a detailed description of the method 
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of solution and presents some comparisons of calculated results with 
exper imental data. 

SYMBOLS 


a 


B 


b 


c 

P 


c 

V 


g»g 1 


H 

h 

j 

M 

P 


r,4> 


b 


00 


00 


speed of sound, a/V^ 
body bluntness parameter 
specific heat at constant pressure 
specific heat at constant volume 
parameters defined by equations (16) 


total enthalpy, H/V^ 

- -2 

static enthalpy, h/V 

00 

indicator, j=0 for 2-D flow, j=1 for axisyrametric flow 
Mach number 

pressure, p/p 

polar coordinates (See fig. 1.) 
gas constant for CF^, 94.475 J/kg-K 

radius of curvature of nose, m 


time, V t/R. 

a> b 

dimensional time, sec 

free-stream velocity, V /V = 1 

J 00 00 


dimensional free-stream velocity, m/sec 
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(V ) 

n s 


v , v A 
r <(> 


v /V 

P 00 


velocity component normal to shock wave, (v ) /V 

ns® 

velocity components in polar coordinate system (see fig. 1), 
and v./V 

<j> ® 


x,y Cartesian coordinates (see fig. 3), x/R b » y/ R b 


Xp location of pole of coordinate system (see fig. 3). X p /R b 

(5 shock wave angle 

s 

6 normal distance between body and shock, m 

Y ratio of specific heats, c /c 

K p v 

Ar local shock layer thickness in radial direction (eq.(12c)) 

An step size in rrdirection 

A£ step size in {(-direction 

e artificial viscosity 

5 angle between tangent to body and line r = Constant 

C angle between tangent to shock and line r = Constant 

s 

n ,n derivative of n with respect to r and 4> (eqs.(12a) and (12d)) 

r <(. 

£,n transformed coordinates defined by equations (9) and (10) 

£ ,£ derivatives of Z, with respect to r and <f> (eqs. (12b) and (12e)) 

r <j> 

p density, p/p 


Superscripts: 

( ) dimensional quantity 

( ) predicted quantity 

Subscripts: 
b body 
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shock 


w wall 

°° free stream 

max maximum 


METHOD 


This section presents a description of the time-dependent, finite- 
difference technique used to solve the inviscid two-dimensional or 
axisymmetric flow over a blunt-nosed body such as that shown in figure 1 . 
The domain of the solution includes the entire subsonic portion of the flow 
field with the downstream outflow boundary located in a region where the 
flow is completely supersonic. 

Flow-Field Equations 

The partial differential equations governing the two-dimensional or 
axisymmetric flow of an inviscid, nonconducting fluid in polar coordinates 
(fig. 1) are 
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denotes the substantial derivative and j=0 for two-dimensional flow and j=1 
for axisymmetric flow. 

Since entropy is constant along streamlines, the following relation can 
be written between pressure and density: 
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( 4 ) 
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which can be used in the continuity equation (1) to replace the density 
derivative to obtain the following equation: 
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Now equations (2), 


(3), and (5) contain derivatives of p, v , 


and v, 
4 > 


only 


and are in a convenient form to solve for these variables. One other 
equation is needed to solve for a second thermodynamic variable. Since our 
interest is only in the steady flow solution, the integrated form of the 
energy equation is used 


H 


h + 



( 6 ) 


Although this will cause the time-varying solution to be inconsistent, it 

will yield consistent and accurate steady-state results. Now equations (2), 

(3). (5), and (6) contain six dependent variables (p, p, a, v , v , and h). 

r $ 

To solve for these variables, two additional equations are needed. These 
can be obtained from the thermodynamic equation of state in the following 
functional form: 


p = p(p,h) 

(7) 

a = a(p,h) 

(8) 


The method of computing these properties for different gases will be 
discussed in the section on thermodynamic properties. 


Computational Domain 


To make integration of this system of partial differential equations easier 
and to help with the application of boundary conditions, the region (fig. 1) 
bounded by the stagnation streamline, the bow shock wave, a line <J> = 4> 

max 

(located in the supersonic flow region), and the body surface is mapped into 
a rectangular domain (fig. 2) using the transformation equations 
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From these transformation equations, a set of transformation operators can 
be defined as 
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which will transform the governing partial differential equations (2), (3) 
and (5) into the following system 
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where 
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Numerical Procedure 


Interior grid points .- At all interior grid points, the governing 
differential equations (eqs. (13) to (15)) are integrated in time, using the 
Brailovskaya difference scheme (ref. 12). This same technique has been used 
by Barnwell (ref. 9) for invisicd flow calculations and Carter (ref. 13) for 
viscous flow calculations. It is an explicit, predictor-corrector technique 
which uses forward time and centered spatial derivatives. 

When it is applied to a general partial differential equation of the 

form 


3f 

at 


+ 



(17a) 


it yields the following result: 


Predictor step: 


f n+1 = f n _ AJL f n ( r n 
i i 2Ax i k i+1 



(17b) 
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Corrector step: 




At Gn *n-1 
2Ax i ' 1 i-1 


f n+1 ) 
i+1 ' 


(17c) 


where n indicates the time step and i indicates the spatial grid point. 

Differencing at the boundaries .- At the shock, body surface, and outflow 
boundaries, three-point, one-sided differences are used for the spatial 
derivatives in a direction normal to the boundary (see fig. 2) in both the 
predictor and corrector steps. Thus, the same order accuracy is maintained 
at the boundaries as at the interior points. Symmetry conditions are 
applied along the stagnation streamline. 

Shock-wave solution .- A method similar to that used by Hamilton and 
Graves (ref. 10) has been used to track the movement of the shock wave 
during the transient portion of the solution. The method can be summarized 
as follows. (Details are given in the appendix.) First, the pressure is 
computed on the downstream of the shock wave by using the two-step 
difference scheme described previously. With the pressure at the shock wave 
P s known, the density is computed for an ideal gas from the equation 
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(18) 


and then the shock velocity is computed from the equation 


v 
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(v) 
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where v g and (v ^ are defined in the sketches in the appendix. Next, the 

movement of the shock wave is computed from the following differential 
equation: 


dr 

s 

dt 
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where 
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. -1,1 dr s. -1 ,S dr s, 

c 3 = tan ( r w ] • tan §r ; 


(21 ) 


by using the following two-step finite difference approximation: 


Predictor step: 


( r b p t + /cos? ) t 

s s s s 


( 22 ) 


Corrector step: 


(r ) t+At “ r o + 0.5At((V /cos? ) t + (V /cos? ) t+At ) 
s s s s s s 


(23) 


For CF^ and equilibrium air chemistry, the Y used in equation (18) is 

replaced by Y . (See section on Thermodynamics.) Thus, for these cases the 

method of tracking the transient shock wave is approximate during the early 
part of the solution, but becomes more exact as the steady state is 
approached. 

Derivatives along the shock, dr /d£, are calculated using a noncentered 

s 

four-point formula of the form 


dr 

( dT } i “ 6AC (r s._ 2 




+ 2r 

s 



(24) 


where i indicates the grid point location. Numerical experimentation has 
shown that this form of the shock derivative produces a much smoother 
converged shock than centered differences. 

Wall boundary conditions .- For inviscid flow, the velocity must be 
tangent to the wall (i.e. the normal component of velocity must be zero). 
This requires that the following equation be satisfied 


(v ) 

r w 


( v ) (tan?. ) 
<p w D 


( 25 ) 


where 
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(26) 


_i i 3p k 

tan (— r— ) 
r b 3 * 


tan 




The wall boundary condition is applied using the following procedure: 

1. The pressure at the wall (p ) is first determined by integrating 

w 

equation ( 13 ). 

2. The "estimated" values for the wall velocity components, 

* * 

(v ) and (v ) , are obtained by integrating equations (I 1 *) 
r w 9 w 

and (15), respectively. 

3. The "final" value for the ^-component of velocity at the wall, 

(v.) , is obtained by combining equation (25) with the 

(p W 

expression for total velocity to obtain: 

(V ) = «vV + (vV)/(1 + tan 2 c )) 1/2 (27) 

<{> w r w w b 


4. The "final" value for the r-component of velocity at the wall, 

(v ) , is then calculated from equation (25). 

5. Using these values of velocity, the wall enthalpy, h , is 

w 

calculated from equation (6). 

This procedure is applied at the wall in both the predictor and corrector 
steps of the integration process. 

Smoothing function .- Numerical experimentation has shown that the 

present solution converges much more rapidly and is less likely to diverge 

for a given set of initial conditions if a small amount of explicit 

numerical smoothing is applied during the integration. In the present 

program, the fourth-order smoother of Barnwell (ref. 9) has been used. The 

smoothing is applied to the dependent variables (p, v , and v) in equations 

r <p 

(13) through (15) during the the corrector step. Small values of the 
smoothing coefficient of the order of 0.02 are typically used in the results 
presented in the present paper and have been found adequate to damp out 
unwanted oscillations without distorting the computed flow field. 

Stability .- The flow-field equations are integrated in the present 
method using the maximun allowable "local" time step at each grid point. 

This has been found to produce almost identical results as solutions 
obtained using the minimum "global" time step, but the speed of convergence 
is "approximately" doubled. To estimate the maximum allowable local time 
step the following approximate CFL relation is used 
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(28) 




a (— + 
^Ar 


1 > 1/2 
rA<p 


where ^ 3 ) m ^ n is the minimum local grid spacing. In actual applications, 

approximately 50 to 80 percent of the time step calculated from the above 
equation is used in the integration at each grid point. 


Thermodynamics 


To compute an actual flow field, density (p) and speed of sound (a) must 
be determined as a function of pressure and enthalpy (p and h) . This is 
expressed in functional form by equations (7) and (8). For an ideal gas, 
these expressions are simply 


P = (fTf) h (29a) 

a 2 = (Y - 1)h (29b) 


Such simple closed form expressions are not possible for CF^. For this 

case, the expressions given by Sutton (ref. 14) have been used to obtain a 
curve fit of 


T = T(h) 


over range of temperatures from 100° to 1000° Kelvin. Thus, once h has been 

computed, T can be determined from the above curve fit, and p can be 
computed from the thermally perfect equation of state: 


p = pRT 


(30) 


where R is the gas constant for CF^ (94.465 J/kg-K). The speed of sound (a) 
can then be determined from the expressions given by Sutton in (ref. 14). 

For air in chemical equilibrium, simple closed form expressions for p 
2 

and a are also impossible to obtain. For this case, Kenneth Sutton has 
developed a set of thermodynamic subroutines based on the approximate 
approach of Hanson (ref. 15) which yield all of the desired thermodynamic 
properties for air at high temperatures as a function of pressure and 
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temperature. To use these subroutines in the present program, an iterative 
scheme was developed to enable all of the desired thermodynamic properties 
to be determined as a function of pressure and enthalpy. 

However, for both CF ^ and equilibrium air, the computational time is 

greatly increased above that required for an ideal gas if the thermodynamic 
routines are called for each time step. Thus a relatively simple procedure 
has been used which requires calling the thermodynamic subroutines only one 
time each 25 to 50 time steps. This procedure requires calling the 
thermodynamic subroutines at the first time step and evaluating the values 
of Y and y a at each grid point which satisfy the equations: 


P = 




(31 ) 


2 

a 


(Y a - 1>h 


(32) 


These values of Y and Y are stored and used in successive time steps to 

pa 

2 

calculate the required values of p and a for each integration step. Since 

2 

the flow-field quanities change slowly, these approximate values of p and a 
are sufficiently accurate to proceed with the solution. At regular 
intervals, the thermodynamic subroutines are called to up date the values of 
Yp and Y^. As the result of several tests, using this procedure requires an 

update to the thermodynamics every 50 iterations. This reduces much of the 
computational overhead so that these types of solutions require only 15 to 
20 percent more time than an ideal gas solution. 


Body Geometry 

The forward portion of the body geometries considered in the present 

paper have a projection in the x,y plane that is described by the equation 
for a general conic section: 


-2 

y 


2R h x - b k 
b b 


-2 

x 


(33) 


where R fe is the radius of curvature at the stagnation point (where x - 0) 

and is a bluntness parameter. The bluntness parameter B^ 

charac ter izes the eccentricity of the conic section. Its significance is 
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0 


better understood if it i.3 noted that B b < 0 generates a hyperbola, - 
generates a parabola, and B^ > o generates an ellipse (with * 1 for the 
special case of a circle). For an ellipse, B b is related to the ratio of 
the major to minor axes (b/a) by the equation 


B fe = (b/a) 2 


(34) 


Some of the body shapes considered in the present paper were blunted cones. 
These shapes can be generated by using equation (33) to describe the nose 
followed by a straight line segment to describe the remaining downstream 
segment of the body. 


RESULTS AND DISCUSSION 


In this section, results of the present method are compared with 
experimental data and the results of other numerical techniques. Two gases 
are considered: air with either ideal gas and equilibrium chemistry, and 

CF^. Flow over a sphere, a paraboloid, and two different ellipsoids is 

considered for the ideal gas cases; flow over a sphere-cone for the 
equilibrium air case; and flow over a sphere for the CF^ case. For all 

cases, Ml points are used around the body and 21 points are used between the 
body and shock. 


Air 


Surface pressure distributions on a sphere in air for Mach numbers of 

2.0, 4.0, and 8.06 and Y » 1 .4 are given in figure 4. To maintain a 

supersonic outflow boundary condition, <t> was set equal to 45 degrees for 

max 

the two higher Mach number cases and 55 degrees for the Mach 2 case. 
Agreement between the present method and the data of references 2, 5, and 16 
is excellent. Shock-layer thicknesses for each Mach number are given in 
figure 5. Again, there is excellent agreement between the present method 
and the numerical and experimental data at the higher Mach numbers. 

However, at a Mach number of 2.0, there is a maximum difference between the 
present method and the experimental data of reference 16 of 6 percent. As 
noted in reference 5, the shock-layer thickness changes rapidly at low Mach 
numbers; and thus, a small error in the test Mach number could be enough to 
cause this difference in results. 

Pressure distributions on a paraboloid for Mach numbers of 3 and 10 and 
Y =* 1 .4 are given in figure 6. The data of the present method are in very 
good agreement with the numerical results of references 2 and 5 for both 
Mach numbers. Shock-layer thicknesses are compared in figure 7, and again, 
agreement is very good. 
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Surface pressure distributions for ellipsoids with b/a =0.5 and 1.5 are 
given in figures 8 and 9 for Y = 1.4 and Mach numbers of 4.0 and 8.06. 
Agreement between the present method and the data of references 5 and 16 is 
excellent. Shock-layer thicknesses for the case where b/a = 1.5 are given 
in figure 10. Agreement is very good for both Mach numbers. 

In figures 11 and 12, surface pressure distributions and shock-layer 
thicknesses are given for flow over a spherically blunted cone, where 
equilibrium air chemistry is used. Comparisons are made with the 
computational data of references 17 and 18. Agreement is generally very 
good except for the shock-layer thickness for values of s/R b > 2.5, where the 

present method produces results that are slightly higher than the other data 
shown. 


The surface pressure distribution on a sphere in CF^ is given in figure 

13. The computer code described in reference 17 was used to generate the 
data for comparison. The agreement is excellent. The computed shock 
standoff distance around the body is shown in figure 14 along with data 
generated by the code of reference 17. Again, the agreement is excellent. 


CONCLUDING REMARKS 


A time-asymptotic method has been used to obtain steady-flow solutions 
for axisymmetric , inviscid flow over several blunt bodies including spheres, 
paraboloids, ellipsoids, and spherically blunted cones. Comparisons with 
experimental data and results of other computational methods have 
demonstrated that accurate solutions can be obtained using this approach. 

The method should prove useful as an analysis tool for comparing with 
experimental data and for making engineering calculations for blunt reentry 
vehicles . 
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APPENDIX 


CALCULATION OF PROPERTIES AT SHOCK WAVE 

A method similar to that presented in reference 10 has been used to 
calculate the shock velocity and other thermodynamic properties at the shock 
wave from the pressure. Consider the shock wave illustrated in the sketch 
(a) which can be expressed as a function of <j>, r - r ($) 

3 3 



The shock-wave angle 6 is given by the equation 

s 


- 2 - ♦ + ^s 


(AD 


where 


- 1,1 3r s, -1 3 v 

«, ■ tan ( r w ) ’ ta " > sr’ 

s s 


( A2) 


and the components of free-stream velocity tangent and normal to the shock 
wave are given, respectively, by the equations 


(V. ) - V cosS 

X, 00 ® s 


(A3) 


( V ) - V sinB 

n ® • s 


(AH) 
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Changes in properties across the shock wave can be related to the 
normal component of free-stream velocity and the shock velocity, as 
illustrated in sketch (b). 



Sketch (b) 


Conditions across the shock can be related by the following normal shock- 
wave equation (since p - 1): 

QD 


(V ) - V 

n ® s 


P„[(V) - VI 

s n s s 


(A5) 
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Combining this equation with equation (A5) and solving for [(V ) - V ] 2 , 

n 00 s 

the following result is obtained 


C(V ) 

n « 


V ] * 



(A9) 


Similarly, equation (A6) becomes 


h 3 ■ h . * it<v. - v 1 c - r ]2 


(A10) 


Equations (A9) and (A10) can be combined with the ideal-gas equation of 
Y - 1 

state p - p( — - — )h to yield the following result: 


h . * j ( P 3 - Pj 
p s < r^r ) * 2 (p s ■ p .> 


(All) 


Using this result, equation (A9) can be solved for the shock velocity V 


P 3 “ P, 2 

V - (V ) - ( — — . — ) 

s n « 1_ 

’ p s 


(A12) 


Using equation (AM), (v ) is obtained 

n s 


(V ) = —[(V ) - V ] + V 

ns p L n ® s s 
s 


(A13) 


Thus the components of velocity downstream of the shock wave are simply 


(v ) = -(V ) cose + V cosp sine 

r s ns s ® s *s 


(AIM) 
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(A15) 


(v ) = (V ) sine + V cosg sine 

<f>s n s s * s s 


For CF^ or equilibrium air chemistry, the Y appearing in equation 


(All) is replaced by 


Y . 
P 


This allows the same shock wave routines to be 


used for these gases as for an ideal gas. Although this approach is 
inconsistent in the transient sense, it becomes more and more exact as the 
steady-state solution is approached. It also reduces the amount of work 
required to obtain converged solutions for these cases. 


18 



nrnnnmrnr'o 
rvhr c*nc*iMOOO 


1. Hayes, Wallace D. and Probstein, Ronald F.: Hypersonic Flow Theory. 

Volume I - Inviscid Flows. Second ed., Academic Press, Inc., 1966. 

2. Van Dyke, Milton D. and Gordon, Hellen D.: Supersonic Flow Past a 

Family of Blunt Axisymmetric Bodies. NASA TR R— 1 , 1959. 

3* Fuller, Franklyn B.s Numerical Solutions for Supersonic Flow of an 

Ideal Gas Around Blunt Two-Dimensional Bodies. NASA TN D-791 , 1961. 

4. Lomax, Harvard and Inouye, Mamoru: Numerical Analysis of Flow 

Properties About Blunt Bodies Moving at Supersonic Speeds in 
Equilibrium Gas. NASA TR R-204, 1964. 

5. Hamilton, H. Harris II: Solution of Axisymmetric and Two-Dimensional 

Inviscid Flow Over Blunt Bodies by the Method of Lines. NASA TP 1154, 
1978. 

6. Moretti, Gino and Abbett, Michael: A Time-Dependent Computational 

Method for Blunt Body Flows. AIAA J., vol. 4, no. 12, Dec. 1966, pp. 
2136-2141 . 

7. Barnwell, Richard W.: A Time-Dependent Method for Calculating 

Supersonic Blunt-Body Flow Fields With Sharp Corners and Embedded 
Shock Waves. NASA TN D-6031 , 1970. 

8. Moretti, Gino and Bleich, Gary: Three-Dimensional Flow Around Blunt 

Bodies. AIAA J., vol. 5, no. 9, Sept. 1967, pp. 1557-1562. 

9. Barnwell, Richard W.: Three-Dimensional Method for Calculating 

Supersonic Angle-of-Attack Flow About Axisymmetric Blunt Bodies With 
Sharp Shoulders and Smooth Nonaxisymmetric Blunt Bodies. NASA TN D- 
6283, 1971. 

10. Hamilton, H. Harris, II and Graves, Randolph A., Jr.: Application of a 

Numerically Generated Orthogonal Coordinate System to the Solution of 
Inviscid Axisymmetric Supersonic Flow Over Blunt Bodies. NASA TP 
1619, 1980. 

11. Weilmuenster , K. James and Hamilton, H. Harris, II: Calculations of 

Inviscid Flow Over Shuttle-Like Vehicles at High Angles of Attack and 
Comparisons With Experimental Data. NASA TP 2103, 1983. 

12. Brailovskya, I. Yu.: A Difference Scheme for Numerical Solution of the 

Two-Dimensional, Nonstationary Navier-Stokes Equations for a 
Compressible Gas. Sov. Phys. - Doklady, vol. 10, no. 2, Aug. 1965, 

pp. 107-110. 

13. Carter, James E. : Numerical Solutions of the Navier-Stokes Equations 

for the Supersonic Laminar Flow Over a Two-Dimensional Compression 
Corner. NASA TR R-385, 1972. 


19 


14. Sutton, Kenneth: Relations for the Thermodynamic and Transport 

Properties in the Testing Environment of the Langley Hypersonic CF., 

Tunnel. NASA TM 83220, 1981. 

15. Hansen, C. Frederick: Approximations for the Thermodynamic and 

Transport Properties of High-Temperature Air. NASA TR R-50, 1959. 

16. Belotserkovsky, 0. M. , ed.: Supersonic Gas Flow Around Blunt Bodies - 

Theoretical and Experimental Investigations. NASA TT F-453, 1967. 

17. Sutton, Kenneth: Characteristics of Coupled Nongray Radiating Flows 

With Ablation Product Effects About Blunt Bodies During Planetary 
Entries. Ph.D. Thesis, North Carolina State Univ. , 1973. (Available 
as NASA TM X-72078.) 

18. Inouye, Mamoru; Rakich, John V.; and Lomax, Harvard: A Description of 

Numerical Methods and Computer Programs for Two-Dimensional and 
Axisymmetric Supersonic Flow Over Blunt-Nosed and Flared Bodies. NASA 
Technical Note, NASA TN D-2970, 1965. 


20 



X 


Shock 



n 



21 



22 


1.2 



Figure 4. - Pressure distribution on a sphere in air. 
ideal gas with y * 1.4. 
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Figure 4. - Continued. 
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(c) ^ = 8.06. 


Figure 4. - Concluded. 
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Figure 6. - Pressure distribution on a paraboloid in air. 
Ideal gas with y * 1.4. 
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(b) M = 10. 
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Figure 6. - Concluded. 
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ShockHayer thicknesses for a paraboloid in air. 
Ideal gas with y =1.4. 
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(a) M =4 
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Pressure distribution on an ellipsoid in air. 
Ideal gas with y =1.4 and b/ a = 0. 5. 
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Figure 8. - Concluded. 
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Figure 9. - Pressure distribution on an ellipsoid in air. 
Ideal gas with y -1.4. and b/a-1.5. 
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Figure 9. - Concluded. 
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Figure 10. - ShockHayer thicknesses for an ellipsoid in air. 
Ideal gas with y - 1.4 and b/a-1.5. 
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Figure 11. - Pressure distribution on a 40° sphere-cone in air 
with equilibrium chemistry. 
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Figure 13. - Pressure distribution on a sphere in CF 4 . 
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